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Preface 


Forward-swept  vinged  aircraft  have  long  been  considered  because 
of  their  aerodynamic  advantages.  However,  until  the  advent  of 
composite  materials,  this  configuration  had  been  structurally  infeasible. 
Now  that  the  forward-swept  wing  is  practical,  the  critical  failure 
mode,  which  past  research  has  found  to  be  divergence,  must  be  explored. 
This  thesis  determines  flutter  speeds,  depending  on  different  aero¬ 
dynamic  models  and  compares  them  to  ascertain  the  practicality  of  more 
advanced  aerodynamics. 

I  would  like  to  thank  Dr.  Robert  Calico,  my  thesis  advisor,  for 
the  time,  help,  and  guidance  that  he  gave  me.  Thanks  are  due  Dr. 

Franklin  Eastep  for  the  aid  given  me  in  the  aerodynamic  modeling  and 
also  Dr.  Wilhelm  Ericksen  and  Mr.  Thomas  Noll  for  their  assistance 
in  the  actual  analysis  and  computation. 

William  L.  Shelton,  Jr. 
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List  of  Symbols  (Ref  l) 


a, 8  Torsional  angular  displacement  of  entire  wing  about  its  elastic 
axis,  in  radians. 

a(s),6(s)  Laplacian  operator  of  torsional  angular  displacement. 

b  Semi-chord  of  wing  including  chord  of  aileron  and  tab,  in  feet. 

h,v  Bending  displacement  of  wing  in  direction  perpendicular  to 
freestream  velocity,  positive  being  down  for  h  and  up  for  w, 
in  feet. 

H(s),W(s)  Laplacian  operator  of  wing  bending  displacement. 

I  Moment  of  inertia  per  unit  span  about  elastic  axis,  in  slug-feet2 
per  foot. 

J=i  V-l 


K  ,Kg  Torsional  spring  constant. 

Kfa,Kw  Bending  spring  constant. 

L  Total  oscillatory  aerodynamic  lift  acting  on  wing  per  unit  span, 
in  pounds  per  foot. 


m  Mass  of  wing. 


M 


Total  oscillatory  aerodynamic  moment  acting  on  wing  about  its 
quarter-chord  point  per  unit  span,  in  pound-feet  per  foot. 


li  Ratio  of  mass  of  wing  to  mass  of  cylinder  of  air  circumscribed 
about  wing  chord,  both  extending  over  equal  lengths.  Note; 


r  Radius  of  gyration  of  wing  about  elastic  axis  as  fraction  of  wing 
semi-chord. 


p  Density  of  air,  in  slugs  per  cubic  foot. 

s  Laplacian  variable  for  characteristic  equations. 

Uj  Critical  divergence  speed,  in  feet  per  second. 

Uf  Critical  flutter  speed,  that  at  which  wing  system  oscillations 
occur  at  constant  amplitude,  in  feet  per  second. 
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U  Forward  speed  of  aircraft  or  freestream  velocity,  in  feet  per 
second. 

Uncoupled  natural  torsional  frequency,  in  radians  per  second. 

Uncoupled  natural  tending  frequency,  in  radians  per  second. 

x  {2(center  of  gravity )/100}-l-a  where  center  of  gravity  is  given 
in  percent  of  total  chord  measured  from  wing  leading  edge.  Note; 


Differentation  with  respect  to  time. 
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Abstract 


Flutter  speeds  for  straight  wings  using  quasi-steady  and  an 
approximation  to  unsteady  aerodynamics  are  determined.  The  same 
procedure  is  then  followed  for  both  straight  and  forward-swept  wings. 
The  equations  of  motion  for  a  rigid-body  aircraft  with  the  proper 
aerodynamics  are  nondimensionalized  and  the  Laplace  transforms  are 
taken.  After  solving  the  coupled  equations  simultaneously,  various 
velocities  are  chosen  and  the  resulting  characteristic  equations  are 
factored.  The  roots  are  then  plotted  and  the  flutter  or  divergence 
speeds  determined. 

It  was  found  that  flutter  speed  increased  and  divergence  speed 
decreased  with  increasing  forward  sweep  as  expected.  Also,  quasi¬ 
steady  aerodynamics  do  not  accurately  predict  flutter  and  the 
approximation  to  unsteady  aerodynamics  does  not  accurately  predict 
divergence.  It  is  necessary  to  combine  these  two  methods  in  order 
to  obtain  a  complete  analysis  of  the  wing's  aeroelastic  instabilities. 


FLUTTER  PREDICTION  OF  FORWARD- SWEPT  WINGS 


BY  ASSUMED  MODES  AND  STRIP  THEORY 


I.  Introduction 

Swept  wings ,  both  forward  and  aft,  have  numerous  advantages. 
Sweeping  the  wing  reduces  the  velocity  component  of  the  freestream 
that  is  normal  to  the  leading  edge,  consequently,  the  swept  wing 
sees  a  reduction  in  drag  and  weaker  shocks  in  supersonic  flight 
(Ref  11:1-3). 

Swept  wings  drew  serious  consideration  back  in  19^7  when 
Collar  studied  their  aeroelastic  phenomena  at  high  speeds.  Diederich 
and  Budiansky  did  the  first  comprehensive  treatment  of  divergence 
in  slender  swept  wings  during  191*8.  They  found  that  the  divergence 
speed  drops  rapidly  as  the  amount  of  forward  sweep  increases  and  that 
wings  with  a  moderate  to  large  amount  of  aft  sweep  cannot  diverge. 

Pros  and  Cons  of  Forward-Swept  Wings 

The  forward-swept  wing  is  not  a  new  concept.  A  four  engine 
Jet  bomber,  the  Junkers  JU-287,  with  15  degrees  forward  sweep,  was 
built  towards  the  end  of  World  War  II.  This  aircraft  actually  flew 
successfully  before  y  ig  damaged  on  the  ground.  More  recently,  the 
German  Hamburger  Flugzeugbau  HFB-300  Hansa  business  Jet  was  designed 
with  a  swept-forvSrd  configuration  to  create  a  more  spacious  fuselage 
interior.  General  Dynamics,  Grumman  and  Rockwell  International  have 
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all  researched  the  swept- forward  wing  and,  in  comparing  it  with  the 
aft-swept  wing,  found  that  it  had  a  lower  stall  speed,  permitting 
higher  manuverahility  at  greater  angles  of  attack  and  improved  slow- 
speed  handling  qualities.  It  exhibited  much  higher  lift  to  drag 
ratios  in  manuvering  flight  as  well  as  lower  trim  drag,  significantly- 
increasing  aircraft  range  in  the  supersonic  regime.  A  wing  combining 
these  characteristics  is  one  having  a  lower  weight  and  therefore 
produces  a  lower  cost  aircraft  for  the  same  high-performance  missions 
then  for  an  aft-swept  wing  made  of  conventional  metallic  materials. 

Until  recently,  the  aerodynamic  advantages  of  forward-swept 
wings  could  not  be  exploited  due  to  its  aeroelastic  instabilities. 

When  a  forward- swept  wing  is  subjected  to  aerodynamic  load,  the  angle 
of  attack  increases.  This  in  turn  increases  the  wing  loading  and 
bending,  which  leads  to  divergence.  To  combat  this  problem  of 
divergence,  the  wing  must  be  stiffened.  However,  this  solution  is 
prohibitive  in  that  the  required  stiffening  produces  a  weight  penalty 
which  negates  any  aeroelastic  advantages.  This  problem  was  studied 
by  Krone  (Ref  6:126-127).  His  report  "clearly  shows  that  the 
detrimental  effect  of  divergence  on  swept- forward  airfoils  can  be 
controlled  by  use  of  advanced  composite  materials"  (Ref  6:126-127). 

Thesis  Objectives 

This  thesis  examines  the  forward-swept  wing,  using  quasi-steady 
and  an  approximation  to  unsteady  aerodynamics,  to  determine  the  wing's 
critical.' flutter  speed.  Tht  ,e  cases  are  compared  and  the  difference 
in  flutter  speed  noted.  The  magnitude  of  this  difference  determines 
whether  both  types  of  aerodynamics  are  used  to  ascertain  the  forward- 
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swept  wing's  critical  flutter  speed.  Once  flutter  or  divergence 
speeds  are  obtained,  the  zeros  of  the  characteristic  equations  are 
traced  in  an  effort  to  gain  some  insight  into  future  flutter 
compensation. 

Solution  Approach 

The  solution  approach  is  a  straight-forward  one.  The  equations 
of  motion  are  written  and  the  proper  aerodynamics  substituted  in. 
Then,  using  strip  theory  and  assumed  modes,  the  resulting  coupled 
lift  and  moment  equations  are  applied  to  the  wing  of  interest.  The 
coupled  equations  are  nondimensionalized ,  then  solved  simultaneously 
through  Laplace  transforms  and  matrix  techniques.  The  coefficients 
of  the  resulting  characteristic  equation  are  functions  of  velocity 
only.  Varying  the  velocity  generates  new  characteristic  equations, 
the  roots  of  which  are  plotted  in  a  root  locus  technique  to  determine 
the  wing's  critical  flutter  speed. 


II.  Theory 


Flutter,  a  dynamic  instability,  is  a  self-excited  vibration 
of  a  body  that  is  moving  in  a  fluid  stream.  Classically,  this  is 
a  mechanism  involving  the  coupling  of  at  least  two  degrees  of  free¬ 
dom  of  the  system.  These  degrees  of  freedom  are  stable  when 
uncoupled  such  that  they  sustain  an  oscillating  motion.  They  are 
critically  coupled  when  there  is  a  phase  difference  in  the 
oscillatory  motion  and  a  time  difference  between  the  aerodynamic 
loads  and  the  system's  resulting  motion.  Flutter  occurs  at  the 
lowest  airspeed  and  corresponding  circular  frequency  that  the 
system,  flying  at  a  prescribed  atmospheric  density  and  temperature, 
exhibits  sustained,  simple  harmonic  motion.  The  forces  necessary 
to  produce  this  motion  are  simply  those  present  when  an  elastic 
structure  is  deflected  from  its  undeformed  state  (Ref  2:527). 

Divergence  differs  from  flutter  in  that  divergence  is  a  zero 
frequency  or  static  aeroelastic  instability  of  an  airfoil  in  torsion 
and  twist.  This  occurs  when  the  torsional  stiffness  of  the  airfoil 
is  exceeded  by  the  aerodynamic  twisting  forces  acting  upon  it  (Ref 
10:298). 

Consequences  of  Instabilities 

These  aeroelastic  phenomena,  when  encountered  unknowingly, 
can  have  drastic  consequences.  Divergence  can  shear  a  wing  .off  an 
aircraft  with  little  or  no  warning.  Flutter  can  result  in  complete 
structural  failure  within  a  matter  of  seconds  or  a  few  oscillations. 
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One  of  the  more  famous  and  spectacular  examples  of  flutter  is 
the  Tacoma  Narrows  Bridge  in  the  state  of  Washington.  The  bridge 
was  constructed  between  two  cliff  faces,  funneling  the  air  between 
them.  When  the  wind  velocity  reached  h2  miles  per  hour,  flutter 
occurred.  This  concrete  and  steel  structure  underwent  displacements 
of  approximately  UO  feet,  behaving  as  though  it  were  rubber. 

Finally,  on  7  November  19^0,  the  bridge  totally  destroyed  itself 
after  undergoing  half  an  hour  of  oscillations  (Ref  7:66). 

Causes  of  Instabilities 

Aeroelastic  phenomena  in  straight  wings  is  dependent  only 
on  the  twist  about  the  elastic  axis.  When  the  wing  is  swept  in 
either  direction,  wing  bending  becomes  very  important  and  greatly 
complicates  matters.  Flax  and  Broadbent  at  the  Royal  Aircraft 
Establishment  indicate  that  for  swept  wings,  the  rigid  body  motions 
are  potentially  dominating.  Flutter  coupling  will  be  between  the 
first  normal  elastic  mode  and  rigid  body  pitching.  In  this  case, 
the  flutter  circular  frequency  will  be  lower  then  the  fundamental 
normal  mode  frequency.  British  investigators  recommend  that,  in 
flutter  analysis  of  swept  wings,  the  first  two  or  three  elastic 
normal  modes  and  two  rigid  body  modes,  pitch  and  translation  for. 
symmetric  modes  and  pitch  and  roll  for  antisymmetric  modes,  be 
used  (Refs  3  and  5)» 

In  a  forward-swept  wing  the  center  of  gravity  is  generally 
aft  of  the  elastic  axis.  This  produces  a  nose-up  torsional  moment 
about  the  elastic  axis  due  to  the  amount  of  lift  produced  by  the 
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airfoil  at  its  aerodynamic  center.  As  the  ving  of  Fig.  1  is  deflected, 
notice  that  line  segment  AB  is  level  while  line  segment  CB  is  canted 
upwards;  point  C  is  higher  then  point  B.  Therefore  segment  BC 
produces  a  positive  increment  of  lift.  This  in  turn  amplifies  the 
torsional  moment  already  existing  about  the  wing's  aerodynamic  center. 
An  unstabilizing  influence,  the  nose-up  twist  shifts  the  center  of 
pressure  outboard.  The  possibility  of  divergence  increases,  but 
the  potential  for  control  reversal  diminishes.  At  speeds  below  what 
is  to  be  called  the  critical  divergence  speed,  these  increments  of 
lift  and  twist  become  smaller  and  smaller  so  that  an  equilibrium 
condition  is  attained.  However,  above  this  critical  speed,  the 
process  is  decidedly  unstable. 

Typical  Sections 

For  typical  multispar,  thick  skin  structures,  it  has 

been  found  that  design  stress  considerations  in  the  vicinity  of  the 
wing  root  results  in  a  root  section  so  stiff  that  the  wing  behaves  as 
a  beam  clamped  perpendicular  to  its  elastic  axis  (Ref  10:3^7). 
Theodorsen  and  Garrick  used  representative  sections  specialized  for 
torsion  and  bending  degrees  of  freedom  (Ref  2:533)*  The  authors 
suggest  that,  for  purposes  of  flutter  prediction,  the  typical  section 
can  be  made  to  fairly  well  represent  a  straight  wing  of  large  span 
by  giving  this  section  the  geometric  and  inertial  properties  of  a 
station  that  is  three-quarters  of  the  way  fom  the  wing's  centerline 
to  its  tip.  This  is  demonstrated  in  Fig.  2. 
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Flutter  Prediction 

Assumed  Modes.  In  predicting  flutter,  the  continuous  aircraft 
ving  is  typically  modeled  using  assumed  modes.  This  procedure 
assumes  a  solution  of  the  boundary-values  associated  vith  a  conservative 
system  of  the  form 

n 

v(x,t)=  T  G  (x)q. (t)  (1) 

i=l  1 

where  G^(x)  is  an  admissible  function,  that  is  one  that  satisfies 
only  the  geometric  boundary  conditions  of  the  differential  equation 
of  the  system,  and  q^{t)  are  the  generalized  coordinates,  to  be 
determined,  which  measure  the  contribution  of  each  of  the  admissible 
QT»  functions  to  the  total  displacement  of  the  wing  surface. 

i  There  are  two  basic  guidelines  to  be  followed  when  selecting 

functions  to  be  used  as  assumed  modes.  They  are; 

1)  Assumed  modes  should  satisfy  the  geometric  boundary 
conditions  and, 

2)  they  should  be  linearly  independent  of  each  other. 

It  is  not  absolutely  necessary  for  the  assumed  modes  to  satisfy 
the  free-edge  natural  boundary  conditions  because  the  process  will 
tend  to  satisfy  them  in  the  final  result.  Linear  independence 
requires  that  it  is  impossible  to  express  one  function  as  the  super¬ 
position  of  the  other  functions. 

Once  assumed  modes  have  been  selected,  the  kinetic  and  potential 
energies  can  be  determined.  It  can  be  shown  that  the  kinetic  energy 

is 
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(2) 


n  n 

KtH  i  z  M.,4.(t)4.(t) 

i*l  J=1  1J  1  1 

and  the  potential  energy  is 

n  n 

V(t)^5  I  r  K  q  (t)q  (t)  (3) 

i=l  J*1  1  1 

where  M^are  the  constant  symmetric  mass  coefficients  and  are  the 
constant  symmetric  stiffness  coefficients,  both  of  which  depend  on 
the  mass  and  stiffness  distributions,  respectively,  and  the  admissible 
functions,  G^(x).  The  potential  energy  generally  contains  spatial 
derivatives  of  G^(x)  that  are  of  orders  as  half  as  large  as  the  order 
of  the  differential  equation  of  the  continuous  system  that  is  being 
considered.  Note  that  the  natural  boundary  conditions  are  of  no  con¬ 
cern  because  they  are  accounted  for  automatically  in  the  kinetic 
and  potential  energies.  Consequently,  G^x)  needs  to  be  an  admissible 
function  only. 

Lagrange's  equations  may  now  be  written  and  are  given  by 


_d 

dt 


9T  I  3T  3  V  „ 


(M 


where  the  q^  represent  the  non-conservative  forces  present.  Using 
Eqs.  2  and  3,  the  equations  of  motion  become 
n  n 

Z  M.q.(t)  +  Z  K  q  (t)  *  0,  r=l,2,...,n 
J«1  rJ  ^  j=l  rj  J 

which,  in  matrix  form,  is 


(5) 


|m|{^(t)}  +  Jk|(q(t))  *  {0} 


(6) 


These  linear  constant  coeficient  equations  can  now  be  solved  for  the 
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q^(t)  and  hence,  along  with  the  assumed  modes  G^(x),  form  w(x,t) 

(Ref  8:266-268). 

Strip  Theory.  Strip  theory  is  used  to  find  the  approximate 
solutions  to  two  dimensional  theoretical  aerodynamic  problems.  It 
assumes  that  the  loads  of  each  spanwise  station  along  the  wing  is 
dependent  only  upon  the  motion  of  that  station  at  any  given  time. 

This  theory  also  assumes  that  there  is  no  spanwise  flow  along  the 
wing. 

The  wing  is  divided  into  "strips"  and  the  aerodynamic  forces 
upon  each  "strip"  are  calculated.  The  aerodynamic  loads  are  based 
on  the  two  dimensional  coefficients  evaluated  at  the  centerline  of 
each  section. 

This  estimation  is  well-suited  for  the  evaluation  of  aerodynamic 
loads  on  a  straight  wing.  When  this  procedure  is  applied  to  a  wing 
with  sweep,  an  approximate  correction,  cosA,  must  be  applied,  where 
A  is  the  quarter-chord  sweep  angle  (Ref  9:17.5-8  -  17.5-9). 

Equations  of  Motion.  When  a  wing  is  represented  by  an  elastic 
surface  that  is  thin  compared  to  its  span  and  chord,  the  equations  of 
motion  can  be  found  through  the  use  of  influence  functions.  The 
wing  is  also  assumed  to  be  rigid  in  the  chordwise  direction  and  is 
then  taken  to  be  straight  and  to  have  no  structural  discontinuities 
so  that  the  existence  of  an  elastic  axis  may  be  assumed.  This  eliminates 
the  elastic  coupling  between  the  integral  equations  of  motion,  leaving 
them  only  inertially  coupled.  Elimination  of  this  inertial  coupling 
can  occur  only  if  the  center  of  gravity  of  each  chordwise  section  lies 
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on  the  axis  (Ref  2:102-105). 

It  is  mor-e' convenient-  to  utilize  the  differential  equations  for 
the  case  described  above.  These  can  be  derived  by  setting  the 
system  to  its  equilibrium  conditions,  then  balancing  its  forces  and 
moments.  Fig.  3  illustrates  the  applied  and  inertial  forces  acting 
on  the  previously  described  straight  wing.  From  this,  the  force  per 
unit  length  is 

Z{y,t )=F  (y,t)+m(y)s(y)e(y,t)-m(y)v(y,t)  (7) 

Zt 

and  the  torque  rate  of  change,  with  respect  to  y,  about  the  elastic 
axis  is 

T,(y,t)*-t(y,t)+s(y){m(y)s(y)6(y,t)-m(y)w(y,t)}+Icg(y)e(y,t)  (8) 

as  seen  from  Fig.  1*  also.  Because  this  is  a  thin  wing,  rotary  inertia 
and  shear  deformation  may  be  neglected.  Therefore,  these  results 


may  be  substituted  into 

{EIv”(y,t)}”=Z(y,t)  (9) 

{GJe’(y,t)),=T,(y,t)  (10) 

Consequently,  the  partial  differential  equations  of  motion  are 

m(y)w(y,t)-Sy(y)e(y,t)+{EIw'’(y,t)}"sFz(y,t)  (ll) 

I  (y)e(y,t)-S v(y)w(y,t)-{GJ6'(y,t)}’*t(y,t)  (12) 

y  * 


•where  is  the  6tatic  mass  moment  per  unit  length  about  the  elastic 
axis  and  I  is  the  mass  moment  of  inertia  about  the  elastic  axis  (Ref 

y 
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Figure  3.  Airplane  wing  with  rigid  chordwise  segments. 
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2:102-106). 

Referring  to  Figs.  5  and  6,  Eqs.  10  and  11  can  be  used  to 
describe  a  beam  under  generalized  loading.  Substituting  the  aerodynamic 
loads  into  the  right  side,  and  taking  El,  GJ,  m,  etc.  to  be  constant, 
produces 


mv(y,t)-S^e(y,t)+EI~trw(y,t)*L(y,t) 

(13) 

Iy6(y,t)-Syw(y,t)-Gj|pr  (y,t)*My(y,t ) 

(HO 

where  L(y,t)  and  M  (y,t )  are  the  aerodynamic  lift  and  moment.  These 
are  the  equations  of  motion  used  in  this  thesis  (Ref  2:5U6). 

Solution  Approach.  Once  the  equations  of  motion  are  obtained, 
several  substitutions  are  necessary  to  put  them  in  a  convienent  form 
to  use.  Notice  that  there  are  fourth  and  second  order  partial 
derivatives  in  the  lift  and  moment  equations,  respectively.  To 
simplify  these  partial  differential  equations,  the  expression  for  a 
beam  in  bending  (Ref  8:209), 

EI-£ppw(y,t)-rou>2v(y,t)*0  (15) 

and,  in  torsion  (Ref  8:19^-195), 

GJ~ye(y ,t )+iyu£  ®(y,t)*o  (l6) 

is  used.  This  yields  the  simplified  the  partial  differential 
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equations  of  motion,  namely 


the  equations  are  nondimensionalized.  Then  the  proper  aerodynamics 
for  each  case  are  inserted  into  the  right  hand  side.  When  assumed 
mode  methods  are  used,  the  bending  mode  is  represented  by 

wty.O.C^)2^)  (20) 

and  the  torsional  mode  by 

e(y,t)=C2sin  ^  ^ei(y,t)  (21) 

vhere  C^and  CgSre  dimensionless  amplitudes  and  (y/b)2  and  sin  ~ 
were  previously  G^(x)  and  w^(t)  and  0^(t)  were  q^(t).  The  lift  and 

moment  equations  are  then  weighted  by  multiplying  them  each  by 


( 


■’»  ■-’•■T  jv.'T;rri> 


* 


5 


^ ;  ■ 

1 

*  © 


I 


I 


C1(^')2  f  22) 

and 

C2Sinfb  (23) 

respectively. 

Strip  theory  is  then  used  to  evaluate  the  aerodynamic  loading 
on  the  right  hand  side  of  the  equations  of  motion.  The  wing  is 
divided  into  five  sections  of  width  0.2b,  where  b  is  the  half-span 
of  the  wing.  Once  this  is  accomplished,  the  loads  on  each  section 
are  integrated  over  the  section’s  width  and  the  terms  summed, 
generating  the  expressions  for  lift  and  moment.  Each  term  on  the  left 
hand  side  of  the  equations  is  also  integrated  across  the  half-span. 

Geometric  and  inertial  values  that  are  constants  for  the  wing, 
straight  or  swept,  are  substituted  into  the  weighted,  integrated 
equations  of  motion.  The  Laplace  transforms  of  these  differential 
equations  are  then  taken  and  put  in  a  matrix  of  the  form 

[A(s)]{xMb}  (24) 

and  solved  simultaneously  by  setting  the  determinant  of  [a(s)]  to  zero. 
The  resulting  polynomial  in  s  has  coefficients  that  are  functions  of 
velocity  only.  Varying  the  velocity  produces  new  characteristic 
equations  which, when  factored,  produce  roots  for  a  root  locus  plot 
similar  to  Fig.  7. 

Because  flutter,  by  definition,  occurs  as  simple  harmonic 
motion,  the  critical  flutter  speed  occurs  as  the  root  locus  crosses 
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the  imaginary  axis.  When  the  real  roots  cross  the  imaginary 
axis,  divergence  occurs  because  this  is  the  zero  frequency 


instability. 
A  through  D, 


Details  of  these  calculations  are  shovn  in  Appendices 
depending  upon  the  case  of  interest. 


III.  Results 


Using  the  methods  described  in  Appendices  A  and  C,  flutter  and 
divergence  speeds  vere  found  for  straight  and  swept  wings  using 
quasi-steady  aerodynamics. 

The  roots  of  the  characteristic  equation,  resulting  from  the 
quasi-steady  analysis  for  a  straight  wing,  are  plotted  in  Fig.  8. 

From  this  root  locus  plot,  the  flutter  speed,  Uf,  and  divergence 
speed,  Uj,  can  be  found.  They  are  the  value  of  velocity  at  which 
the  high  frequency  branch  crosses  the  imaginary  axis  for  flutter  and 
the  low  frequency  branch  crosses  the  imaginary  axis  for  divergence. 
Therefore,  from  the  root  locus,  it  can  be  seen  that  Uf  is  62.5  ft /sec 
and  Uj  is  217.8  ft/sec.  Notice  the  inordinately  low  flutter  speed. 
Because  this  is  an  invalid  value,  as  explained  later,  root  locus 
plots  will  not  by  done  for  the  cases  of  15  and  30  degrees  of  forward 
sweep.  The  constant  term  of  the  characteristic  equation  is  solved 
for  velocity  in  the  above  two  cases  and  the  divergence  speeds  are 
found  to  be  157.7  ft/sec  for  -15  degrees  of  sweep  and  109.2  ft/sec 
for  -30  degrees  of  sweep. 

The  exceptionally  low  flutter  speeds  are  a  result  of  the  quasi¬ 
steady  analysis  being  inadequate  for  use  at  higher  reduced  frequencies. 
Guidelines  for  use  of  the  appropriate  aerodynamics  with  their  corres¬ 
ponding  frequency  ranges  can  be  ascertained  through  the  use  of  reduced 
frequency. 

The  reduced  frequency  is 
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Figure  b.  Hoot  locus  of  flutter  analysis  of  a  straight  wing  with 
quasi-steady  aerodynamics. 
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Axes  have  units  of  rad/sec 
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,  o»b 

k=r 

To  find  the  flutter  reduced  frequency,  use  the  flutter  speed  and 
frequency.  For  quasi-steady  aerodynamics ,  k  must  be  less  than  0.2. 
This  reduced  frequency  restriction  is  rather  strict  and  is  rarely 
met  by  typical  systems  at  subsonic  speeds.  As  seen  in  Table  1,  none 
of  the  quasi-steady  aerodynamic  cases  are  near  this  value,  therefore 
these  are  not  valid  flutter  analyses.  Because  divergence  is  a 
zero  frequency  instability,  quasi-steady  aerodynamics  predict  this 
phenomena  quite  well.  Therefore  the  divergence  speeds  are  valid 
values  and  are  plotted  in  Fig.  9  as  a  function  of  sweep  angle,  A. 

The  methods  of  Appendices  B  and  D  utilize  an  unsteady  aero¬ 
dynamic  model  to  find  the  flutter  and  divergence  speeds  for  a 
straight  wing  and  ones  that  are  swept  forward  15  and  30  degrees. 

The  unsteady  aerodynamics  utilized  an  approximation  to  the  comp¬ 
licated  Theodorsen  functions.  This  was  Jones's  approximation  and 
is  valid  only  at  higher  reduced  frequencies.  The  roots  of  the 
characteristic  equations  for  each  of  these  three  cases  are  plotted 
in  Figs.  10,  11  and  12,  respectively.  From  Fig.  10,  is  shown 
to  be  165.3  ft/sec.  For  15  and  30  degrees  of  forward  sweep,  Uf 
is  l8l.2  and  201.1  ft/sec,  respectively.  Note  that  the  divergence 
speed  cannot  be  determined.  In  the  unsteady  case,  the  reduced 
frequency  range  changes,  as  seen  in  Table  1,  and  is  above  the  lower 
limit  of  0.2.  Consequently,  the  flutter  speeds  predicted  by  the 
unsteady  aerodynamics  approximation  are  good  values.  The  variance 
of  the  flutter  speed,  as  predicted  by  unsteady  aerodynamics,  is 
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gure  10*  hoot  locus  of  flutter  analysis  of  a  straight  wing  with 
unsteady  aerodynamics. 
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Figure!!;.  Root  locus  of  flutter  analysis  of  &  wing  rvept  forward 
15  degrees  using  unsteady  aerodynamics. 
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Table  1.  Ranges  of  Reduced  Frequency 


Reduced  frequency,  k 

Sweep  angle,  A 

0° 

O 

ITS 

H 

1 

-30° 

Quasi-steady 

aerodynamics 

1.11 

0.95 

0.90 

Unsteady 

aerodynamics 

0.29 

0.25 

0.21 

shown  in  Fig.  15*  The  same  trend  in  flutter  in  flutter  speed  is 
displayed  here  as  it  was  in  the  quasi-steady  aerodynamic  cases. 
Un fortunately,  the  approximation  to  the  unsteady  aerodynamics  do 
not  do  a  good  job  at  predicting  the  critical  divergence  speed. 
Again  this  is  due  to  the  discrepancies  in  the  reduced  frequency 
range . 


IV.  Conclusions 


Flutter  analysis,  using  assumed  modes  and  strip  theory,  is  a 
straight  forward  process.  To  get  valid  results,  however,  the 
proper  aerodynamics  must  he  used  for  each  case. 

Quasi-steady  and  the  approximation  to  the  unsteady  aerodynamics 
are  relatively  easy  to  use  and  do  not  require  the  use  of  a  computer. 

Once  the  unsteady  aerodynamics  are  sweep-corrected  though,  it  is 
impractical  to  perform  the  necessary  algebraic  manipulations  by 
hand.  Even  though  the  use  of  the  MAC  Symbolic  Manipulation  System, 
see  Appendix  D,  reduces  this  tedious  job  significantly,  it  is 
still  highly  inefficient.  Therefore,  unless  it  is  unavoidable, 
these  expressions  should  not  be  used. 

The  two  aerodynamic  cases  examined  here  are  at  the  opposite 
ends  of  the  frequency  spectrum.  This  leaves  a  gap  which  potentially  may 
be  filled  by  tht  Pade  approximation.  This  approximation  represents 
the  transformed  aerodynamic  forces  as  a  ratio  of  polynomials  that 
are  a  function  of  the  Laplace  transform  variable  s,  as  did  the 
Jones's  approximation.  The  coefficients  of  these  polynomials  are 
found  using  a  two  step  least  squares  process.  The  generalized 
aerodynamic  forces  are  obtained  from  the  Pade  approximation  of 
sinusoidal  forces  using  doublet-lattice  aerodynamics  (Ref  : 885-886). 

The  aeroelastic  instability  of  flutter  occurs  when  at  least 
two  degrees  of  freedom  of  the  system  coalesce.  Typically,  for  aft- 
swept  wings,  this  is  between  the  first  bending  and  first  torsion 
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modes.  When  looking  at  forward- swept  wings  the  short  period  pitching 
of  the  aircraft  combines  with  its  low  frequency  mode,  producing 
flutter.  The  Pade  approximation  should  be  compared  with  the  quasi¬ 
steady  analysis  and  that  of  the  unsteady  approximation  to  find  the 
theory  that  yields  the  most  accurate  results. 

This  thesis  has  laid  the  groundwork  for  this  comparison  and 
an  analysis,  such  as  the  one  described  above  for  a  swept  wing, 
would  be  a  logical  follow-on. 
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Appendix  A 

Flutter  analysis  of  a  straight  wing  using  quasi-steady 
aerodynamics  is  a  straight-forward  process.  The  equations  of  motion 
are  (Ref  2:5^6) 


4 


mtt(y,t)-S  9(y,t)+EI^-vw(y,t)=L(y,t) 

“  3y 

2 

(Al) 

I  6(y,t)-S  V(y,t)-GJ— y9(y,t)*M (y,t) 

°  a  3y  y 

(A 2) 

The  equations  for  a  beam  in  bending  and  in  torsion,  respectively,  are 
(Ref  8:209,  19U-195) 

94  2 

El — qv(y,t)-niD  v(y,t)=0  (A3) 

3y 

32  2  , 

GJ — z$(y,t)+I  u>  9(y,t)=0  [Ah) 

3y  00 

Substituting  the  expressions  for  quasi-steady  aerodynamic  lift  and 
moment  (Ref  2:279), 

L(y,t  )*2TtpUb{U9(y,t)-w(y,t)}  (A5) 

My(y,t)eL(y,t)b(a+^)  (a6) 

and  Eqs.  A3  and  Ah  into  Eqs.  A1  and  A2,  respectively,  results  in 
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mV-S  0+mo),  v=2TTpUb(U0-v) 
cx  n 


2  2 

I  e-S  V+I  OJ  e=2TrpUb  (a+^}(ue-v) 
a  a  a  a 


To  make  use  of  the  method  of  assumed  modes,  each  equation  of  motion 
must  he  weighted  by  its  generalized  coordinates.  The  functions 
C^^-)  and  Cgsin  ^  were  chosen  to  weight  the  lift  and  moment  equations, 
respectively,  because  the  bending  of  the  wing  is  closely  approximated 
by  a  parabolic  and  the  twisting  by  a  sinusoidal  function. 

For  similar  reasons,  the  assumed  mode  shapes  of  bending  and 
torsion  are  taken  to  be 


w(y,t)sC1(^-)2w1(t) 

8(y,t)=C2sin|^ei(t) 


(A10) 


When  Eqs.  A5  and  A6  are  weighted  and  mode  shapes  substituted  in,  Eqs. 
All  and  A12  are  the  results. 

Pim(b)\cic2s«(t)2sinib5i',cim“h(b)\- 


C.2irpb(£)  U(UG-w) 

±  D 


(All) 


si,,2ibei’ 


C22rpUb2(a+H)sin|^-(  U6-w) 


(A12) 


Unploying  strip  theory  requires  that  the  right  hand  sides  of 
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i 

Eqs.  All  and  A12  "be  divided  into  "strips"  and  integrated  over  that 
strip.  Strips  of  width  0.2b  were  chosen.  The  left  hand  side  must 
also  be  integrated,  but  this  is  done  only  across  the  semi-chord,  b, 
and  not  divided  into  strips.  Once  this  has  been  accomplished, 

Eqs.  All  and  A12  are  nondimensionalized  using  the  quantities. 


S  =mbx 
a  a 


m 


I  =mb2r2 
a  a 


(A13) 


V 


( 


Once  this  has  been  accomplished,  the  weighted,  integrated 
equations  are 


O.2J1-  £l(0.295)x  e.40.2agjl- 
d  Cn  a  1  no 


1.257  1  U  v, .1.851  C„  1  ,U*A 
v  V  b  b 1*T-  cl  v  (b)ei 


(All*) 


6.-  ~1  ( 0 . 589 )~ar,l+uj2 0 . ■ 
!  C2  rab  a  1 


U 

V 


7^ 


{0.1*94(^)261-^1(0.295)^  ^1} 


(A15) 


These  equations  must  be  solved  simultaneously  and  to  do  so, 
their  Laplace  transforms  must  be  taken.  It  is  assumed  that  the  initial 
conditions  are  zero  60  consequently,  the  transformed  equations  are 
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{0.2s’*±^pIA£s  ♦  0.2«J)  fils) 


^2(0.295  xos 


2+  aj5i  W)2}6  (s)«o 

TT  V  D  1 


(A16) 


-^1(0.589  f?s 
c2  o 


1.18  (aj^)  U 
it  r  2  b 


s}  il(s) 


+{s2+u)2_  ii22i  Is^(^)2)e  (s)-o  (AIT) 

aw  r*  b  1 

For  the  wing  being  analysed, 

u>  =10  rad/sec 
n 

id  =25  rad/sec 
a 

x  =0.1 
a 

r  =0.5 
a 

a=-0.2 
b=3  ft 
y=20 

Using  these  constants  in  Eqs.  Al6  and  AIT , 

(s2+0.033Us+100)2l-£2(0.lkTs2+O.Ol6U2  >6^0  (Al8) 

-§l(0.236s2-0.02fcUs)2l4(s2+625-0. 013U2)6  -0  (A19) 

C2 

which  can  easily  be  put  into  matrix  form,  the  determinant  of  which 
will  yield  the  desired  characteristic  equation: 
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s4+0.038Us3+(751.06-0.018U2)s2 
+21. 56lUs+(61*T1+6. 71-1. 367U2  )=0  ( A20) 

This  equation  has  coefficients  that  are  functions  of  velocity  only. 
Varying  the  velocity  produces  nev  characteristic  equations  and  the 
roots  of  these  are  plotted  in  a  root  locus. 
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Flutter  analysis  of  a  straight  wing  utilizing  unsteady  aero¬ 
dynamics  follows  the  same  basic  procedure  as  that  for  a  straight  wing 
with  quasi-steady  aerodynamics,  the  basic  difference  being  the  more 
complex  aerodynamic  expression.  The  same  wing  equations  of  motion, 
generalized  coordinates,  and  mode  shapes  from  Appendix  A  are  used. 

The  new  aerodynamic  expressions  are 


L(y,t)=rrpb2{Ue(y,t)-tt(y,t)-bae(y,t)} 

+2irpUbC(k){U0(y,t  )-w(y,t  )+b(5g-a)0  (y,t )}  (Bl) 

My(y,t)=-*pb2{baV(y,t)+Ub(3s-a)e(y,t)+b2  (-g+a2  )e(y,t)} 

+2irpUb2(a+^)C(k){Ue(y,t)-w(y,t)+b(^-a)e(y,t)}  (B2) 


(Ref  2:272)  where  C(k)  is  a  Theodorsen  function  that  is  dependent  on 
the  reduced  frequency,  k.  Therefore,  substituting  Eqs.  Bl  and  B2 
into  Eqs.  A1  and  A 2  from  Appendix  A  results  in 


*  -  X  5  ♦  -  £  -.5)  *  £  2  cUJ#  -  £  ♦<*-.)«! 

ba  hbpbb  y  b  bb 


r2  -x  ^  +  r2u29*  -  ^  (*s-a)0  +(i 

aob  aa  y  b  b  o 


+a2  )0) 


,(a+1l)  U  wU-  w 


Because  Theodorsen  functions  are,  in  general,  very  complicated,  Jones's 
approximation 


1«0 


(B5) 


ff1 


* 


.  5s2+  0.81<3Us  +  O.OlUU2 
~9s2+  1.037Us  +  O.OlMJ2 


(Ref  2)  will  be  utilized  instead. 

Following  the  same  procedure  outlined  in  Appendix  A,  the 
weighted,  integrated  equations  of  motion  are 

VC  **  \T  C  V 

0.2^  -  ±2  0.0295x  6.+0.2ui2  ^1=  ^2  0.295  -  ^0,-0. 2^1 
DC.  a  1  h  b  C,  p  b  1  b 


0.295  J  6  +  §2  0.598  i($C(s)  -  —  ~C(s)Jl 

1  Jj  D  pDD 


+  &  0.598  £  (%-a)C(s)e, 
L1  D 


(B  6) 


°-5rI®r  |  0-?9^ab1W'5l'»“»61-|  °-295  *•  |l 

b  P  1  p  1  b  p  1 

-71  0.591-^^-r  c(8)  ~i  +  g  c(s)  e  (bt) 

C„  p  b  b  p  b  1 


Transforming  Eqs.  B6  and  BT  and  using  the  constant  values  of  Appendix 
A,  yields. 


{0.l*s2+0.0067UC(s)s+20}£l 

D 

-72{0. 032s2+(0. 005+0. 0069C(s))Us-0.0033U2C(s)}6*0  (b8) 

C1  1 

-7l{0.032s2+(0. 0058-0. 0035C(s)s}Jl 
C2  15 

+{0. 129s2+{ 0. 0058-0. 0035C( s ) )Us 

+(78. 125-0.  oorri^cUWe^o 


(B9) 


and  solving  them  simultaneously  produces  Eq.  BIO  which  is  the 
characteristic  equation  of  the  sy;tem: 

s6+0 . 156Us5+( T3U . 2-0. 0031U2 )s4+( 98.75-0. 00132U2 )Us3 
+  ( 60385+2. 62U2-0. 000625U4 ) s  +(6953-0.  O89U2 )Us 
+(92. 07U2-0. 007U4 )=0  (BIO) 

Again,  Eq.  BIO  has  coefficients  that  are  functions  of  velocity 
only.  The  results  of  varying  the  velocity  are  plotted  in  a  root 


locus. 


Appendix  C 


Flutter  analysis  of  a  swept  wing  does  not  differ  greatly  from 
that  of  a  straight  one.  The  major  difference  occurs  in  the  definition 
of  the  variables  that  are  dependent  upon  the  spanwise  axis.  The 
variables  in  question  are  the  semi-chord,  b,  and  the  mass  ratio,  y . 
They  are  defined  as 


b— 

cosA 


itpUb5 


where  A  is  the  sweep  angle,  positive  being  aft.  The  redefinition  of 
the  semi-chord,  b,  also  affects  the  mode  shapes: 


v*C.(^)2w. 

b 


0=C?sin  —  ^  0, 
2  b 


The  quasi-steady  aerodynamic  expressions  for  lift  and  moment 


L(y,t )*Tpb2{U0cosA+UasinA-UbaTsinA} 
+2irpbUcosA{-w+U0cosA+UosinA 
♦b(*s-a)(0+UxsinA)} 


M  (y,t )=-irpb3{U(i5-a)ecosA+3sU2TCOSAsinA 

y 

-UaasinA+b(^+a2 )UtsinA) 

+2irpb2UcosA(35+a){-w+U6cosA 

+UasinA+b(?g-a)(e+UTsinA)}  (C6) 

.  -3v  .  8? 

where  o=  —  and  t=  —  . 

3y  8y 

Eqs.  C5  and  C6  are  the  unsteady  aerodynamic  expressions  that 
are  simplified  by  the  quasi-steady  assumptions,  namely  small  disturbances 
thereby  eliminating  the  accleration  terms  and  C(k)  is  one  (Ref  2:398). 

Because  divergence  is  :a  zero  frequency  aeroelastic  instability, 
only  the  constant  term  of  the  characteristic  equation,  which  has  been 
greatly  complicated  by  the  sweep  correction,  is  of  importance. 

This  term  of  the  weighted,  integrated,  transformed  equation  is  set  to  0: 

U2(0.0007sin2Acos2A-0.99sin22A+0.Usin2A-0.3cosA) 

+U(0.0007sin2A) 

+(-0. 0002sin22A-0 . lsinA-2 . 19sin2A+15625 )»0  ( C7 ) 

Selecting  A  to  be  either  -15  or  -30  degrees  and  solving  for  U  will 
give  the  divergence  speed  for  that  case. 


This  calculation  was  greatly  complicated  by  the  addition  of 
the  sweep- corrected  unsteady  aerodynamic  terms.  The  new  expressions 
for  lift  and  moment,  respectively,  are 


L(  y , t ) ='  pb2  {  W+U6  c o s A+Uo s  inA-ba6 -baUr s inA  } 
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where 
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3y 
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3y 
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y 
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(Ref  2:398). 

Following  the  procedure  of  Appendix  B,  the  nor.dimer.sionalized 
equations  for  lift  and  moment  are 
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mb  pb 
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The  equations  are  weighted  and  integrated  using  the  same 
constants  and  equations  of  motion  as  before,  then  taking  their 
Laplace  transforms,  resulting  in 
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vhere  C(k)  is  the  Theordorsen  function  that  is  approximated  by 
Eq.  B5. 

Solving  these  complicated  expressions  by  hand  would  be  very 
difficult,  if  not  impossible.  Therefore,  this  was  accomplished  by 
using  the  MAC  Symbolic  Manipulation  Syatem  (MACSYMA)  that  was 
developed  by  the  Massachusetts  Institute  of  Technology  Math  Lab. 

This  system  is  an  interpeted  computer  language  for  algebraic 
manipulations. 

When  Eqs.  D6  and  D7  are  solved  using  MACSYMA,  the  expression 
is  quite  long  and  complicated  and  will  not  be  shown  here.  This 
equation  is  a  function  of  sweep  argle  and  velocity,  which,  when  the 
velocity  is  varied  for  angles  of  15  and  30  degrees  of  sweep,  produces 
characteristic  equations.  The  roots  of  these  are  plotted  on  a 


root  locus. 
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